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Abstract 

A method is introduced to calculate the UV-divergent parts at one-loop level in dimensional 
regularization. The method is based on the recursion, and the basic integrals are just the scaleless 
integrals after the recursive reduction, which involve no other momentum scales except the loop 
momentum itself. The method can be easily implemented in any symbolic computer language, and 
an implementation in Mathematica is ready to use. 
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PROGRAM SUMMARY 

Title of program: $UVPart 
Programming language: Mathematica 

Available from: http : //power . itp . ac . cn/~f engf eng/uvpart/ 

Computer(s) for which the program has been designed: Any computer where the 
Mathematica is running. 

External routines/libraries used: FeynCalc, FeynArts 

Keywords: UV-Divergences, One-Loop Corrections, Dimensional Regularization 
CPC Library Classification: 11.1 

Nature of problem: To get the UV-divergent part of any one-loop expression. 

Method of solution: $UVPart is an Mathematica package where the recursive method 
has been implemented. 

Running time: In general it is below a second. 
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LONG WRITE-UP 



I. INTRODUCTION 



One has to deal with an integration over the loop momentum at next-to-leading 
order, which results to ultraviolet (UV) and infrared (IR) divergencies. Dimensional 
regularization[H, 0] is needed in order to produce meaningful results. The general one- loop 
amplitude can be written as of a linear combination of known scalar integrals [3[ - boxes, 
triangles, bubbles and tadpoles — multiplied by coefficients that are rational functions of 
the external momenta and polarization vectors, plus a remainder which is also a rational 
function of the latter. 

There are many automatic tools available to achieve the general one-loop amplitude, like 
FeynCalcPJ and FormCalc^, which are based on the traditional Passarino- Veltman 
0] reduction of Feynman graphs, which can be generated automat ically ( Feyn Art s fiol. 
or QGRAf(12|). In order toproduce numerical results, tensor coefficients functions are 
calculated using LoopToOLsjl]. For a detailed review, please see Refs. fl3l. 

In the last few years, several groups have been working on the problem of constructing 
efficient and automatized methods for the computation of one-loop corrections for multi- 
particle processes. Many different interesting techniques have been proposed: these contain 
numerical and semi-numerical methods [loT- 18l|. as well as analytic approaches (l9 - 22 1 that 



make use of unitarity cuts to build NLO amplitudes by 
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luing on-shell tree amplitudes |23l. 
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For a recent review of existing methods, see Refs. [25|, 
Generally, it will be much easier to calculate the UV-divergent parts of the one-loop 
amplitude than the one-loop amplitude itself, and there are little work on this specific area 
since one usually need not to calculate this part separately. However if what we are concerned 
is the renormalization, we need to calculate the UV-divergent parts only, and moreover there 



are also some cases in which we have to calculate the UV- divergence separately, for example 



most packages like Fire[27| and Reduze[28|, which implement the integration by parts 



(IBP) relations [29|, treat the scaleless integrals as zero, i.e. 
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however it is well known for the logarithmically divergent scaleless integrals that 
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so Eq. ([T]) does not distinguish the UV- and IR- divergence. This will be fine if we consider 
the amplitude as a whole, since the UV-divergent parts will be canceled by the counter- 
terms, and the left divergence will only be IR-divergent. But there is no way to know the 
IR-divergence from a specific Feynman diagram, which is very important when one considers 
the factorization, where one tries to identity the source of IR- divergence and factorize them 
out. 

A method that allows the extraction of the UV-divergent part of an arbitrary 1-loop 
tensor N-point coefficient was presented in Ref. 30(. We want to introduce another simple 



method to calculate the UV-divergent part in the dimensional regularization. The method is 
based on the recursion, and basic integrals are just the scaleless integrals after the recursive 
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reduction, which involve no other momentum scales except the loop momentum itself. Since 
the computation in this method just involves algebraic rational operations, so it can be 
easily implemented in any symbolic computer language. 

The paper is organized as follows: We introduce the definitions and notations in Sec. HTl 
then describe the calculations in Sec. IIII[ and in Sec. IIVI we give an implementation with 
Mathematica, and use this method in a specific process e + e~ — > J/ip + f] c , and finally 
comes the summary. 



II. DEFINITIONS AND NOTATIONS 

The general expression associated with UV-divergent parts at one-loop reads 

_ (2n^ D f Q 

m 2 J N?N^---NT [ } 

with the denominators 

Mi = (k + qi ) 2 - m 2 + ie, (4) 

where ie denotes the infinitesimally small positive imaginary part, /i is the reorganization 
scale, D is the non-integer dimension of space-time defined as D = 4 — 2e, qi are linear 
combination of external momenta pi, and the numerator Q is the polynomial of k 2 and k-pi. 

The following identity about the scaleless integrals is well known in the calculations with 
dimensional regularization: 

(2^)"- J />,, ,,. 2r JO (n/-2) 



where only the logarithmically divergent scaleless integral contributes the UV-divergence. 
According to Lorentz covariance and oddness of the scaleless integrals, we have 

/ d o k w . . . k , m {kT = ( (m is odd) 

in 2 J I Crn,n x g { ^ 2 Mm| (m is even) ' v ' 

i.e. when m is odd, the scaleless tensor integrals are 0, and when m is even, they are 
proportional to g^ lf12 ' "^™1 with coefficient C mi „, where we use the same notations as Refs. jl, 
30| for g{wf ,2 '"' 1 "'} ) which is the symmetrization of g^ 1 ^ 2 g^ 4 ■ ■ ■ g^m-iiMn w jth re spect to the 



lorentz index /i 2 , ■ ■ • , /i m , for example 

^{/ilM2^3M4} _ glHH2 g^3^4 _|_ glUlLZ gillie _|_ ^Ml^^^MS ^ 

can be ready achieved by mult 
SW2SW4 ■ ••9n m -mm and contracting the indexes, 



ffl 

The coefficient C m , n can be ready achieved by multiplying Eq. ((6]) with the — metric tensors 



c - I (2?r/i)4 D f d D k(k 2 Y + ^ 

Lm ' n ~ D(D + 2) ■■■(D + m- 2) ztt 2 7 1 > 
(n + 2 ^ -2) 
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where we have used 

W W 4 • ■ ■ =D(D + 2) ••■ (D + m - 2) 

where m is even, and more relations can be found in Ref. (30| . 



(9) 



III. DESCRIPTION OF THE CALCULATIONS 

We define a function Power to get the asymptotic scaling of the polynomial of k 2 and 
k ■ pi in the ultraviolet region, for example: 



PowerfA; 2 ] = 2, Po-wer[k ■ p^ = 1, Power[m 2 ] = 0, Power[fc 2 + k ■ pi 
and we can extend this function to rational expression of k 2 and k ■ Pi. 

"AT 



m 



Power 



V 



Power [N] - Power [£>], 



where M and D are some polynomials of k 2 and k ■ Pi , for example 

k-pi 



Power 
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So if Power of the integrand in X is less than -4, then there will be no UV-divergent part in 
X. i.e. 



iV 
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uv 



For example 
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(k 2 — ml) ((k + pi) 2 — mf) 2 

Now we are going to describe the calculations of UV-divergent parts of X. For general 
integrand of X, we can write 
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with Q' = Q (A/i — k 2 ), and it is clear that 
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i.e. either the power of propagators: rii or the Power of the integrand decreases by at least 
1, if we apply this replacement once again in the last result, we have 
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with Q" = Q! [M\ — k 2 ), and we can see the power of propagators: rii or the Power of the 
integrand decreases further by at least 1. 

So we can apply this replacement again and again until one of the following cases happens: 

• All the power of Mi, becomes zero, i.e. only one type of propagator: k 2 is left. 

• The corresponding Xuv = 0, i.e. 



Power 



Q 



M^M^-^-M, 



n N 
N 



N 



Power[Q] - < -4. 



So after the recursive reduction, only one type of integration will be left: 



X 



d D k 



Q 



and it is ready to read the result according to Eqs. ([5]) and ([6]). 
We can take the following 3-point tensor integral as an example: 

(k-px) 3 
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(18) 



(19) 



where each step — > means we take a replacement, and all the expressions framed with box 
have been dropped during the recursive expansion since the Power is less than — 4 and will 
not contribute UV-divergent part. 

Now it is ready to get the UV-divergent parts from the last line in Eq. (Tl9|) . using Eqs. ([5]) 
and (El) 



A-D 



d°k 



(k-Pi) 
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and we can check that it agrees with Refs. 

la lag. 

The important feature of this method is that, there is only one preliminary integral, i.e. 
Eq. ([6]), which involves no other scales like external momenta pi or mass m^. 

Another advantage is that it can be easily implemented in any symbolic computer lan- 
guage, like Mathematica, Reduce, Form, etc. I will give an explicit implementation 
with Mathematica in the next section. 



IV. IMPLEMENTATION WITH MATHEMATICA 



A. Typical Examples 



An implementation in Mathematica is already available, note that the FeynCalc 
package has been used to deal with the ScalarProduct, but it not required for the imple- 
mentation. 

The UV-divergent parts can be retrieved with the function: $UVPart, which is defined 

as: 

<f°k exp 



$UVPart [exp , k] : = 



(2tt/x) 4 



(21) 



uv 



where exp can be any expression at one-loop level, while k is the loop momentum. 
We can give some simple examples: 

den=SPD [k] (SPD [k+ql] -ml~ 2) (SPD [k+q2] -m2"2)//FCI//ScalarProductExpand; 
num=SPD[k, pl]~3//FCI; 
$UVPart [num/den , k] //Simplify 



The output of the code above reads: 



Out[ ] := 



6pl 2 (pl -ql +pl -q2) 
D(D + 2)co 



(22) 



where u is just £uv which we have used to represent the UV divergence in the Mathematic 
code, Eq. (122 j) gives the same result as Eq. (12 Op after setting the dimension D to 4. 



den = (SPD[k]-mO~2) (SPD [k] -ml~2) "2 (SPD [k] -m2~2)//FCI//ScalarProductExpand; 
num = SPD[k,p] "8//FCI; 
$UVPart [num/den , k] 

The output reads: 

105 (m0 4 



0ut[ ] 



setting D to 4 we get 



(2ml 2 + m2 2 ) mO 2 + 3ml 4 + m2 4 + 2ml 2 m2 2 ) (p 



2\4 
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128 
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(24) 
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B. UV-Divergences in e + e — > J/ip + i] c at One-Loop Level 

We will apply the method to a specific process e + e~ — > J/ip + r] c , where we use the 
FeynArts to generate Feynman diagrams, and FeynCalc to handle the Dirac Trace. We 
will take a triangle Feynman diagram in Fig. [TJ as a concrete example in this section. 



■f WWIAAM 




j/i> 



FIG. 1: An triangle feynman diagram in e + e — > 7* — > J/ip + rj c 

First we use the FeynArts to generate the amplitude for the diagram, then make the 
following replacement with the help of spinor projectors j3lj : 



v{p 3 )u{p 3 



1 



-(#, - m c )e£(f 3 + 2^3)^3 + m c ) 



4V2E 3 (E 3 + m c , 
v(p 4 )u(p 4 ) ->■ -(^ 4 - m c ) 75(^4 + 2£ , 4 )(j4 4 + m c ) 



4^2(^4 



1 

W 



(25) 



where the subscripts 3 and 4 are used to label J/ip and 7/ c respectively. 

After performing the DiracTrace on the fermion chains, we get the amplitude for this 
diagram as follows: 



Amp 



6y/W c m 5 c {s - A)s 2 k 2 (k 2 + k- p A ) (2sm 2 + k 2 + k ■ p 3 + 2k ■ p 4 ) 
x {-2Ds 2 m A c + \2s 2 m\ + %Dsm\ - 56sm* + 32m\ + 8Dk 2 m 2 c 



-2Dsk 2 m\ 



Ask 2 m\ 



ry ry ry ry ry 

16k m c + ADk ■ p 3 m c + D sk ■ p 3 m c 



—12Dsk ■ p 3 m 2 c + 36sk ■ p 3 m 2 — 16k ■ p 3 m 2 c + ADk ■ p^m 2 c 
—D 2 sk ■ p^m 2 c + lODsk ■ p^m 2 c — 32sk ■ p^m 2 c — 2Dk ■ p\ 

+Ak ■ p\ + 2Dk ■ p 3 k ■ p 4 — Ak ■ p 3 k ■ p 4 



(26) 



where p 3 and p 4 are the momenta of J/ip and rj c respectively, and s = (p 3 + p 4 ) 2 /(Am 2 ), m c 
is the quark mass, and k is the loop momentum, and to get the UV-divergent part of the 
amplitude, we just use 



$UVPart [Amp,k] 
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The output reads: 



Out[ ] 



iC F (D - 2fegie 



2 P(1 4 7^p3p4 



(27) 



where the Lorentz index 7 refers to /x, and ip to e* s in Eq. ( J26l) . After setting .D to 4, we get 



(2vr/i) 



4-£> 



Z7r 



d D k Amp 



uv 



(28) 



We can apply this method to each diagram to get the corresponding UV-divergent part 
of the amplitude, and to check the validity of the our result. We have compared the UV- 
divergent part produced with our code with Ref. 32| diagram by diagram, and the both 
agree with each other for all diagrams. 



V. SUMMARY 

A pretty simple method is introduced to calculate the UV-divergent parts at one-loop level 
with dimensional regularization. It is found that there is only one preliminary integral which 
involves no other scale like external momenta p,i or mass m, after the recursive reduction. 
The method can be easily implemented in any symbolic computer language, An explicit 
implementation with Mathematica is also present. 
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